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ABSTRACT 

A fraction of the extragalactic near-infrared (near-IR) background light involves redshifted photons 
from the ultraviolet (UV) emission from galaxies present during reionization at redshifts above 6. The 
absolute intensity and the anisotropics of the near-IR background provide an observational probe of 
the first-light galaxies and their spatial distribution. We estimate the extragalactic background light 
intensity during reionization by accounting for the stellar and nebular emission from first-light galax- 
ies. We require the UV photon density from these galaxies to generate a reionization history that is 
consistent with the optical depth to electron scattering from cosmic microwave background measure- 
ments. We also require the bright-end luminosity function of galaxies in our models to reproduce the 
measured Lyman drop-out luminosity functions at redshifts of 6 to 8. The absolute intensity is about 
0.1 to 0.4 nW m~^ sr~^ at the peak of its spectrum at ~ 1.1 /im. We also discuss the anisotropy 
power spectrum of the near-IR background using a halo model to describe the galaxy distribution. 
We compare our predictions for the anisotropy power spectrum to existing measurements from deep 
near-IR imaging data from Spitzer/IRAC, Hubble /NICMOS, and AKARI. The predicted rms fluctua- 
tions at tens of arcminute angular scales are roughly an order of magnitude smaller than the existing 
measurements. While strong arguments have been made that the measured fluctuations do not have 
an origin involving faint low-redshift galaxies, we find that measurements in the literature are also 
incompatible with galaxies present during the era of reionization. The measured near-IR background 
anisotropics remain unexplained with an unknown origin. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium — large scale struc- 
ture of universe 



1. INTRODUCTION 

The optical and UV radiation from sources present 
during reionization is expected to leave a signature in the 
extragalactic background light (EBL) at near-IR wave- 
lengths (e.g., Santos et al. 2002; Salvaterra & Ferrara 
2003; Cooray & Yoshida 2004; Fernandez & Komatsu 
2006; Raue 2009). Such radiation is not expected to be 
present in the background light at UV and optical wave- 
lengths due to the redshifted Lyman limit. The exact 
intensity from first-light galaxies present during reion- 
ization is currently unknown. The first predictions sug- 
gested an intensity as high as 10 to 30 nW m~^ sr~^ 
(Santos et al. 2002; Salvaterra & Ferrara 2003). These 
estimates were partly motivated by the need to explain 
the difference between DIRBE EBL measurements (e.g., 
Cambresy et al. 2001) and the integrated galaxy light 
(IGL) from deep galaxy counts (Madau & Pozzetti 2000; 
Totani et al. 2001). 

These predictions with high backgrounds were ques- 
tioned by Madau & Silk (2005) based on existing limits 
related to metal content at high redshifts and the X- 
ray background produced by stellar end-products such as 
black holes. They suggest an intensity less than about 
2.5 nW m~^ sr~^ in the J-band from a galaxy popula- 
tion made up of population HI stars during reionization 
(Madau & Silk 2005). With a combination of Popula- 
tion II stars and changes to the lifetime of stars, Fer- 
nandez & Komatsu (2006) argued that the background 
could be as high as 4 to 8 nW m~^ sr"'^. Even in such 



a scenario a simple estimate of the UV photon density 
at z > 6 shows that there are roughly an order of mag- 
nitude higher number of H-ionizing photons per baryon 
during reionization than necessary to explain the reion- 
ization history. Since one does not expect more than a 
few H-ionizing photons per baryon during reionization, 
a first-order estimate suggests that the background in- 
tensity cannot be larger than a few tenths nW m~^ sr~^ 
between 1 and 2 /xm. 

Unfortunately a direct search for the integrated inten- 
sity of galaxies present during reionization based on ab- 
solute background measurements has been problematic 
due to the confusion with the Zodiacal foreground. At 
1 AU Zodiacal light is two to three orders of magnitude 
brighter than the ^ 10 nW m~^ sr~^ intensity produced 
by extragalactic sources. While challenging, techniques 
have been devised to estimate the Zodiacal dust column 
density based on the line strengths of Fraunhofer lines 
seen in the dust-scattered Solar spectrum (e.g., Bernstein 
& Dyson 2003). Instead of the absolute background, in 
Cooray et al. (2004; also Kashlinsky et al. 2004), it was 
proposed that the galaxies present during reionization 
can be studied with anisotropics of the near-IR back- 
ground. The anisotropy studies have the potential to 
probe deeper than the absolute experiments and could 
study a galaxy population present during reionization 
that leads to an intensity well below 0.1 nW m^^ sr~^ 
(Cooray et al. 2004). 

This suggestion has motivated experimental measure- 
ments on the near-IR anisotropy power spectrum with 
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data from Spitzer/IRAC, HST/mCMOS, and AKARI. 
After a deep removal of point sources, Kashlinsky et 
al. (2005, 2007, 2012) claimed a detection of first-light 
galaxy fluctuations at z > 8. The detected signal is an 
excess of clustering power above shot-noise on the largest 
angular scales. A similar suggestion was also made by an 
AKARI gronp (Matsumoto et al. 2011), but an analysis 
of the HST/mCMOS Ultra Deep Field led to an opposite 
conclusion that the sources contributing to the near-IR 
excess fluctuations are at z < 8 (Thompson et al. 2007). 
Due to the limited areas of existing deep surveys near- 
IR background anisotropy measurements are limited to 
angular scales less than about one degree. The limited 
field of view is especially a problem for existing NIC- 
MOS UDF measurements (Thompson et al. 2007), where 
the fluctuations are limited to angular scales less than 
5'. Separately, a joint analysis of IRAC and HST/ ACS 
data in the same GOODS fields as studied by Kashlin- 
sky et al. (2007) led to the suggestion that up to 50% of 
the excess fluctuations at 3.6 /im could come from faint 
dwarf galaxies at z ^ 2 (Cooray et al. 2007; Chary et 
al. 2008). Through detailed models combined with more 
recent measurements of faint galaxy clustering, Helga- 
son et al. (2012) has lowered this low-redshift contribu- 
tion to 3.6 ^m intensity fluctuations to be at most 20%. 
The rest of the anisotropies continue to be interpreted as 
originating from first-light galaxies during reionization 
(Kashhnsky et al. 2012). 

While there are still uncertainties on the exact inten- 
sity and the amplitude of intensity fluctuations in experi- 
mental measurements, the situation is no different on the 
theory side. The first estimates on the anisotropy power 
spectrum made use of linear theory clustering (Cooray 
et al. 2004). Fernandez et al. (2010) used numerical 
simulations of reionization to predict the expected power 
spectrum during reionization. Their power spectrum has 
a shape in the form of a power- law with Ci oc Z^/^ be- 
tween ten arcminute to arcsecond angular scales. Fernan- 
dez et al. (2010) suggested that the power-law behavior 
arises from significant non-linear biasing of dark matter 
halos at high redshifts. Due to the limited box sizes of 
existing reionization simulations of the order 100 to 140 
Mpc on the side, numerical studies are limited to angular 
scales of 30 arcminutes and below at z > 6. 

With the availability of WFC3 on the Hubble Space 
Telescope, dedicated IR background experiments (e.g., 
GIBER sounding-rocket experiment; Bock et al. 2006), 
and plans for a future space-based absolute intensity 
measurement (ZEBRAQ; Cooray et al. 2009) there IS now 
a clear need to revisit theoretical predictions on both the 
absolute intensity and the anisotropy power spectrum 
from galaxies present during reionization. While cur- 
rent multi-band Spitzer and AKARI measurements do 
not overlap in the same fields, the combination of IRAC 
and WFC3 data on some of the same well-studied fields 
on the sky (e.g., fields covered by the CANDELS survey; 
Grogin et al. 2011; Koekemoer et al. 2011) will soon 
allow the spectral energy distribution (SED) of intensity 
fluctuations be studied uniformly. Separately CIBER is 
conducting spectral imaging absolute measurements be- 
tween 0.8 and 1.6 /xm in wide 4 deg.^ fields instanta- 
neously using multiple sounding rocket fiights (Zemcov 

^ http://zebra.caltech.edu 



et al. 2012). The combination of IRAC and CIBER 
is capable of extending anisotropy measurements out to 
angular scales of more than a degree from optical to 4.5 
fim. 

In this work we establish both the mean intensity and 
the anisotropy power spectrum of galaxies present during 
reionization. We update Cooray et al. (2004) by taking 
into account recent developments in the study of reion- 
ization, and by introducing a halo model to calculate the 
non-linear clustering of the IR background intensity. The 
stellar and nebular emission from first light galaxies fol- 
low the calculations presented in Fernandez & Komatsu 
(2006), but we specifically require that the UV photon 
background produced by the galaxy population present 
from z ~ 6 to 30 is consistent with the optical depth 
to electron scattering as measured by the WMAP po- 
larization data with a value of 0.088 ±0.014 (Komatsu 
et al. 2011). We account for the current uncertainty in 
the optical depth by introducing variations to the fidu- 
cial model so that the optical depth to electron scatter- 
ing varies between 0.07 and 0.1. Our models are also 
designed to reproduce the bright-end galaxy luminosity 
functions (LFs) in deep IIST/WFC3 surveys at z > 6 
involving the Lyman-dropout galaxy samples. This nor- 
malization at the bright-end of galaxy luminosities puts 
strong constraints on the intensity. 

This paper is organized as follows. In Section 2 we out- 
line our model for the reionization galaxies including stel- 
lar nebulae and the IGM emission. Section 3 presents the 
calculation related to luminosity mass density of these 
galaxies. In Section 4, 5 and 6 we outline the background 
intensity and spatial anisotropy power spectrum calcula- 
tions, respectively. In Section 7 we discuss our results 
related to the intensity and angular power spectrum and 
present a comparison to existing measurements. We con- 
clude with a summary in Section 8. We assume the flat 
ACDM model with = 0.27, = 0.046, Ha = 0.73 
h = 0.71, as = 0.81 and = 0.96 ([Komatsu et al.ll201lD 
throughout the paper. 

2. EMISSION FROM STARS AND THE INTERGALACTIC 
MEDIUM 

We first describe the emission from stars in first-light 
galaxies present during reionization. Following Fernan- 
dez & Komatsu (2006), we consider two stellar popula- 
tions in this calculation. The first, referred to as Pop II 
stars, are metal-poor stars with metallicity Z = 1/50^0, 
and the second. Pop III stars, are metal-free stars with 
Z = 0. 

To describe the stellar initial mass function (IMF) we 
make use of two descriptions. For Pop II stars, we adopt 
the IMF given by Salpeter (1995) 

f{AU) cx AC'-'', (1) 

with mass range from 3 to 150 Mq. For Pop III stars we 
use the IMF obtained by Larson (1999), which takes the 
form as 

/(^^*) = AC^(l + ^) (2) 

where = 250 Mq, and the mass range is from 3 to 
500 Mq. 

we utilize the fitting results from iLeieune fc Schaereii 
pool.) and .Schaerer, (,2002,) to calculate other stellar 
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parameters, such as the intrinsic bolometric luminosity 
the effective temperature T;^{M^), the main- 
sequence lifetime r*, and the time-averaged hydrogen 
photoionization rate (5jjj(M*). The fitting forms of these 
parameters are different for Pop II and Pop III stars. For 
Pop II stars, they are given as 

logio(L!j°V^0) = 0-138 + 4.28a; - 0.653x^ 
logio(7T^7K) = 3.92 + 0.704a; - 0.138a;^ 
logio(T,/yr) =9.59 - 2.79a; + 0.63a;^ 
logio(QHi/s"^) = 27.80 + 30.68a; - 14.80a;2 + 2.50a;^ 
where x = logio(Af*/Mo) and for Pop III stars, they are 

logio(i^°V^0) =0.4568 + 3.897a; 0.5297a;2, 
Iogio(r5/K) = 3.639 + 1.501a; - 0.5561a;2 + 0.07005x^ 
logio(T,/yr) =9.785- 3. 759x -I- 1.413x2 - 0.186a;3, 



logio((3Hi/s ^) 



39.29 8.55a; 5-9 Mq 

43.61 + 4.90a; - 0.83a;2 9-500 Mq 



From these expressions the stellar radius R^:{M^,) is 



47ri?2(M, 



(3) 



where a — 5.67 x lO^'^ erg s^^ cm^^ K^'' is the Stefan- 
Boltzmann constant. The stellar radius is useful for the 
calculation related to the stellar emission spectrum (see 
Section 3.1). 

The ionization volume in the nebulae surrounding the 
stars (Stromgren sphere) can be derived if assuming ion- 
ization equilibrium where the ionization rate equals re- 
combination rate 



'HII 



^rcc ' 



(4) 



where a^'^ is the hydrogen case B recombination coef- 
ficient which depends on the gas temperature Tgas (as- 
sumed to be ~ 3 X 10^ K), and we will discuss it in detail 
in the next Section. Here, n"'''' and ng|j are the local 
number density of electron and HII in the stellar nebulae 



where we assume n 



neb 



^ncb 
'HII 



lO'' 



For the IGM, the hydrogen density is lower than that of 
the stellar nebulae, so we no longer assume the ionization 
equilibrium. We estimate the ionization volu me by a 
redshift-dependent form as (jSantos et al. l[2002l ) 



IGM, 



mi 



(M,) 



(5) 



where nii{z) — 1.905 x 10 ^(1 + z)'^ cm ^ is the mean 
hydrogen number density for Qi, = 0.046 assumed in this 
work (ShuU et al. 2011). 

These quantities discussed here would now be used 
to estimate the luminosity mass density, near-IR back- 
ground intensity SED, and the anisotropy power spec- 
trum. We make use of emission from the stellar nebulae 
and the IGM for both Pop II and Pop III stars in galaxies 
present during reionization. 



3. LUMINOSITY MASS DENSITY OF THE SOURCES 

In this Section, we calculate the luminosity per stellar 
mass at frequency v, i.e. luminosity mass density, for sev- 
eral sources that contribute to the infrared background, 
such as the direct emission from the stars, Lyman-a line, 
and free-free, free-bound and two photon processes. The 
luminosity mass density takes the central role in our es- 
timation of the near-IR background intensity spectrum. 

3.1. Stellar spectrum 

For simplicity, we assume the stellar spectrum is a 
Planckian truncated at hi' = 13.6 eV. Thus, the stellar 
luminosity at frequency v can be expressed as 







for hiy < 13.6 eV 
for hi/ > 13.6 eV 



where S"* — 47ri?^ is the surface area of the star, i?* is 
the stellar radius, Af* is the stellar mass and B^iT) is 
the Planck spectrum 



bat: 



icff\ 



(6) 



Note that for simplicity we have ignored the absorption 
lines of the Lyman-alpha (Lya) series here. The absorp- 
tion is not strong enough to affect the shape of the spec- 
trum and hence it is not expected to affect our results es- 
pecially in the infrared wavelengths ( Santos et al. |[2002l : 
iFernandez fc Komatsul[2006i) . Our predictions related to 
the optical background fluctuations may be somewhat 
overestimated. Also, the emission with hv > 13.6 eV 
cannot be approximated by a black-body spectrum, thus 
we use the fitting formulae for time-averaged photoion- 
ization rate, Qhij Pop H ^'^id Pop III stars to calculate 
the emission at higher energies. 

3.2. Lyman alpha emission 
The luminosity of Lya emission at a frequency v is 



/l^Lya(eLya + ^Ly" )'^('^Lya " 1^, z)V{M.), (7) 



where z^Lya is the frequency of Lya photons, and V{M^) 
is the emission volume that can be estimated by Eq. Q 
and Eq. ([5]) for the stellar nebulae and the IGM respec- 
tively. Here CLyc is the Lya recombination emission rate 
per cm'^, which is given by 

^rcc jrec ^ ^ ^ rcc 

^Lya = Jhya^enmiaB ) 

where n,, is the electron number density, tihii is the HII 
number density, a^'^ is the hydrogen case B recombi- 
nation coefficient, and /£y'^„ i^ the fraction of the Lya 
photons produced in the case B recombination. This 
fraction can be estima ted through the fitting formula 
(jCantalupo et al.ll2008D 



flTJT) = 0.686 - 0.1061ogio(T4) - O.OOgT^ 



-0.44 



(8) 



where T4 = T/IO^'K, which is accurate to 0.1% for 
100 < T < lO^K. Note that this fraction is actually 
not sensitive to the temperature, so it can be treated 
as a constant ~ 0.68 in the most cases. The hydrogen 
case B recombination coefficient a^ we use here is from 
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Fig. 1. — The luminosity mass density vs. the rest-frame wavelength A for the Pop II and Pop III stars at 2 = 10. The total from 
the stellar nebula and the IGM are shown respectively in each figure (thick lines). Also, to specify the contribution from different sources, 
we plot the of stellar, Lyman-a, free-free, free-bound and two-photon emission for the stellar nebula case (thin lines). We take /* = 0.03 
and /esc = 0.5 for both Pop II and Pop III cases. We find the of the stellar nebula is much larger than that of the IGM for both Pop II 
and Pop III stars. 



iHummeH (|1994D . which is fitted by iSeager et~aLl ()1999f ) 

as 



where 



a'^^T) = 10 



4.309, b 



-13 



1 



[cmVs], 



(9) 



-0.6166, 



0.6703 and d 



0.5300. We assume a gas temperature of Tg, 
K in our calculation to obtain ctg'^. 



3 X 10^ 



The 



by 



^coll 
^Lya 



is the coUisional emission rate per cm given 



coll 



(10) 



where nui is the neutral hydrogen number density and 
C£yQ, is the effective coUisional excitation coefficient. It 
has the form (jCantalupo et al.ll2008l ) 



^Lya 



Cl,2p + Cl, 



3s 



3d ■ 



(11) 



Here we take into account the excitation up to energy 
level n = 3 to produce Ly-a photons. The higher level 
emission can be neglected given the high temperature we 
consider for this calculation. The excitation coUisional 
rate Ci,u, in cm'^ per second, can be written as 

8.629 X 10-6 



C, 



91 



T 



-11, 



,6-^-/'=^ [cmVs], (12) 



where E\^^ is the energy difference between lower level / 
and higher level m, gi is the statistic weight for level I, 
and ji,u(T) is the effective col lision strength calculate d 
using the fitting formulae from iGiovanardi et al.l ()1987f ). 

The 0(i^LyQ — y^z) is the Lyg line profile, and we use 
the result from iSantos et al. I (I2002D wh ere they fitted 
the simulated Lya line profile of iLoeb fc"Rvbicki (1999| ) 
for a homogeneous and expanding IGM: 



{z)dv ^exp[— i/, (z)/(ii'] 



if di^ > 
if di^ < 



where dv = 



and 



,{z) = 1.5 X 10"Hz 



0.019 



h 

OJ 



[l + zf 
E{z) 



Here E(z) = H{z)/Hq, and we assume the fiat ACDM 
model and take E{z) = ^yQ^(l + z)^ + fl^. 

3.3. Free-free and free-bound emission 

For free-free and free-bound emission we again fol- 
low the same approach as Fernandez & Komatsu (2006). 
Following their derivation, the continuum luminosity of 
these two processes at frequency v has the same form 
with 

L^^fb = 4nj«^^V{M,), (13) 

where j^'^ is the specific emission coefficient for free- free 
and free-bound emission UDopita fc Sutherlandl[200^ 



jf'fi^ = 5.44x10" 



39 f 



-hv/kT 



T 



ff,fb 



[erg /cm^ /s/Hz /sr] . 

(14) 

Here Up is the proton number density, T is the gas tem- 
perature, and g^g^ is effective Gaunt factor for free-free 
and free-bound emission, which takes the form as 



ff,fb 
ffcff 



9« 



-fffb(n) 



free- free 
free-bound 



(15) 



where gg ~ 1.1 is the thermal averaged Gaunt factor of 
free- free emission and gfh{n) ~ 1.05 is the free-bound 
emission Gaunt factor for a different energy level n. 
These values have an accuracy of 10% (H arzas fc Latteii 
I1961D . In above Xn = Ry/{kTgn'^), where Ry/kTg is 
aro und 10 for the parameter s pace we are interested 
in (|Fernandez fc Komatsul I2006D . The energy level n 
is determined by the emission photon frequency v. If 
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cRy/n'"^ < V < cRy/{n' — 1)^, and then n — n' where 
Ry = 1.1 X 10^ m^^ is the Rydberg constant. Note that 
the n here starts at n = 2, since the photons from n = 1 
can be easily absorbed by other neutral hydrogen atoms 
and be ionized instantly. 

3.4. Two-photon emission 

For the two-photon process we also follow the approach 
of Fernandez & Komatsu (2006) and write the luminosity 
as 

^2ph ^ 2h!Lp^^/^^^^^^^^v{M,), (16) 

where e2pii = /2ph?ie'T-HiiQ!B'^ two-photon emis- 

sion rate per cm^, and /2ph — (1 ~ fhya)- ^he 
P{v I vi,ya)dv I vi^ya IS the normalized probability of gen- 
erating one photon in the range dv/vi^ya from per two- 
photon de cay. We use the fitting formula derived in 
[Ferna ndez fc Komatsul ()2006D 

P{y) ^ 1.307 - 2.627(y - 0.5)^ 

+ 2.563(j/ - 0.5)'' - 51.69(?/ - 0.5)^ 

where j7 = dv/uhva. which is a good fit to the data 
given in iBrown fc Mathews! ([TOTOf ). 

3.5. Luminosity mass density and total emission 

Then following Fernandez & Komatsu (2006) and Fer- 
nandez et al. (2010), we can derive the mean luminosity 
mass density for each emission source by integrating over 
the IMF for the Pop II or Pop 111 stars 



/ dMJ{A'QL^{M, 
(M,) 



(17) 



where the ranges of the integral are from 3 to 150 Mq 
with the IMF given by Salpeter (1955) for Pop II stars, 
and from 3 to 500 Mq with the IMF in Larson (1998) 
for Pop III stars. Here, the average mass (M,) is given 
as 

(18) 



Jo 

where /(Af,) is the normalized IMF with 
dMJ{A'Q = 1. 



Note that this expression is only valid when the main 
sequence lifetime is larger than the star formation time 
scale. Otherwise, it should be evaluated by 



/ dM,f{M,)L^{M,)T,{Ah) 
tsY{z){M,) 



(19) 



Here isF is the star formation time scale which is given 

by 

' d\x\pf,{z) ' 



tsFiz) 



dt 



(20) 



where p[z) is the stellar mass density at z, which is 
related with the comoving star formation rate density 
(SFRD) as -0(2;) = dp^{z)/dt. We use the halo mass 
function to calculate ip{z) and p*(2;) and the details are 
described in the next Section. We note that the tsF is 



important for the esti mation of the 1^ as discussed in 
[Fernandez et al.l (|2010D . 

Finally we obtain the total luminosity mass density 
from the stellar nebulae 

= 11 + {I- /csc)(?J;''" + C + C + iT) (21) 

and the same from the IGM 

z^^^ = /csc(?j;''"+c+e+e'), (22) 

where f^sc is the escape fraction of the ionization photons 
that propagate into the IGM from the nebulae surround- 
ing the stars. 

In Fig. (TJ we show the total luminosity mass density 
as a function of the rest-frame wavelength A from the 
stellar nebula and the IGM for Pop II and Pop HI stars 
respectively. The contributions from the different sources 
we consider are also shown for the stellar nebula case. 
Here we set /, = 0.03 and /esc = 0.5 for both Pop II and 
Pop HI cases. 

With these parameters we find that the stellar spec- 
trum is dominant for Pop II stars while the "background" 
spectrum, such as Ly-a and free-bound, are comparable 
with or even larger than the stellar spectrum for Pop HI 
stars. Also, as can be seen, the lu from the IGM is much 
lower than that from the stellar nebula for both Pop II 
and Pop HI cases, and the total lu from the Pop HI stars 
is similar to that from the Pop II stars. These results are 
already discussed in Fernandez & Komatsu (2006). 

4. REIONIZATION HISTORY AND UV LUMINOSITY 
DENSITY 

To test if the reionization history associated with our 
stellar model is consistent with that of the cur rent obser- 
vatio ns, such as WMAP 7-year results (Komatsu et al.l 
1201 Ih . we need to calculate the hydrogen reionization 
fraction XYai{z) as a function of redshift. Following 
Madau et al. (1998), xy{ii{z) can be estimated as 



dx 



HIT 



fcsc'4'{z)q{z) 



dt 



xmi 



(23) 



nu{z) 

where Zj is the redshift of the beginning of the reion- 
ization epoch (we take Zi — 30), ^{z) is the comov- 
ing star formation rate density (SFRD), and the func- 
tion q{z) is defined as q{z) = (Qhi/(-^^*)) (■^*)- Here 
(r*) is the average stellar lifetime which is given by 
(r*) = dM^ T^{M^)f{A'U), and iroc is the volume 
averaged recombination time, which can be written as 



[CHIl(^)aB'^H(^)(l + l'/4X)]- 



(24) 



where Chii = ("■hii)/("-hii)^ is the clumping factor of 
ionized hydrogen. Here we adopt the simulation result 
from Trac & Cen (2007). X = 0.75 and Y = 0.25 are 
the mass fractions of hydrogen and helium, respectively. 
Note that we have already considered the escape fraction 
/esc, so the Chii here is the clumping factor excluding the 
halos with star formation. 

For the comovi ng SFRD we c onsider the ongoing star 
formation model ([Santos et al. 1 12002'): 



nu d 



1"°° dn 
^i^)^f*i^^J dMM—(M,z), 



dM 



(25) 



where is the star formation efficiency which denotes 
the fraction of baryons converted to stars, dn/dM is the 
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T = 0.077 
T = 0.090 
T = 0.099 

I L 




halo mass function (jSheth fc TormerJ 119990 . and Mmin 
is the threshold mass for a dark matter halo to form a 
galaxy during reionization. This minimum mass is taken 
to be a free parameter and is varied to fit a combination 
of the WMAP 7-year optical depth and the galaxy LF, 
as we discuss in Section 6. 

We also need a stellar population evolution model to 
describe q{z) with the relative fraction of the Pop II and 
Pop III stars at different redshifts. In principle there 
should be a cutoff at some redshift for Pop III stars as 
they are not expected to form at low redshifts once the 
gas is polluted by metals. We assume this cutoff is not 
lower than z = 6 when the universe is fully ionized. We 
use the error function to denote the population fraction 
as 



Fig. 2. — Left: The hydrogen reionization fraction xhii as a function of redshift z for three cases of the /, = 0.02, 0.03 and 0.04. The 
sohd line with /« = 0.03 indicates the optical depth t = 0.090 which is close to the result of WMAP 7-year which gives t = 0.088 it 0.014. 
The dashed and dotted lines are derived with /, = 0. 02 a nd 0.04 which denote t = 0.077 and 0.099, respectively. Right: The dependence 
of the optical depth on the minimum redshift in Eq. I I28II . The blue solid, dashed and dotted lines are for the three stellar models with 
/* = 0.03, 0.02 and 0.04. The result of WMAP 7-year is also shown in red solid line with yellow la error region for comparison. 

is the electron number density of the Universe at red- 
shift z, and we assume the helium is singly ionized for 
z > 4 (77 = 1) and doubly ioniz ed for z < 4 (77 = 2) 
(|Kuhlen fc Faucher-GiguCTglMl ). 

In Fig. [21 we plot the hydrogen reionization fraction 
a^Hii vs. z for three optical depth r with three stellar 
models. The blue solid line denotes the r = 0.090 with 
/* = 0.03 which is close to th e result of WMAP 7- year 
data with r = 0.088 ± 0.014 (jKomatsu et all 120 111) . In 
this case reionization ends around Ze — 8 which is consis- 
tent with the current studies, and we find Pop III stars 
can ionize the Universe to ^ 60% by z = 10 and Pop II 
stars is responsible for the rest of the reionization over 
an interval Az ~ 2. 

For /, = 0.02 and 0.04, we find r = 0.077, shown 
with a dashed line and r = 0.099, shown with a dotted 
line, respectively. These two reionization histories are 
such that Zc ~ 6.5 and 9, respectively. For case with r = 
0.077, the Pop III stars ionize ~ 40% of the Universe and 
the Pop II stars are needed over the interval Az ~ 3.5 to 
complete reionization. For the third case with r = 0.099, 
Pop III stars ionize 80% of the Universe at z = 10 with 
Pop II stars completing the rest 20% in an interval of 
Az ~ 1 

In Fig. [21 we also show the dependence of the optical 
depth on the minimum redshift Zmin in Eq- (|28l) . The 
blue solid, dashed and dotted lines are for the three stel- 
lar models with r = 0.090, 0.077 and 0.099 for Zmi„ = 0. 
The WMAP 7-year resuh with r = 0.088 ± 0.014 is also 
shown with a red solid line and a yellow la error region. 
We find the slope of the curves is steeper for Zmin > 10 
and flatter for z,„i„ < 10 which is caused by the Pop 
III to Pop II transition around 2 = 10 in our model. 
Note that this transition is arbitrarily chosen. We varied 
the transition redshift and also cases where PopII and 
PopIII stars are mixed in with different fractions at dif- 
ferent redshifts. In all these cases we found results that 



1 + erf 



zt - 10 

crp 



(26) 



where crp = 0.5 is the population transition width. Then 
the term q{z) in Eq. (j23p can be expressed by 



q{z) = /p 



-^PopIII 



^PopIII 



>+(l-/p) 



T^Popll 



PopIlN 



(27) 

As our default model we assume that the Pop III stars 
are mainly dominant for z > 10 while Pop II stars are 
for z < 10 with zt = 10. When we present our results 
we also show results for three additional values of the 
transition from zt = 10 to 30. 

We estimate the optical depth to electron scattering 
with the reionization fraction Xe(z) as 



dz 



ne(z)o-T 



H{z) 1 + z 



(28) 



where ut = 6.65 x 10 m^ is the Thompson scatter- 
ing cross-section, and ndz) = xmi{z}nii{z){l + jijY/AX) 
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are generally consistent with each other. Thus, the three 
choices related to the reionization history that we show 
here to keep this presentation simple are not biased with 
respect to the final result related to the IR background 
intensity that we are trying to estimate in this paper. 



CO 
J3 



O 




Fig. 3. — The total number of ionizing photons per baryon to 
maintain the ionized IGM between the redshift of the end of the 
reionization and z. The blue dotted, solid and dashed lines are 
derived from the model with /« = 0.02, 0.03 and 0.04, respectively. 

We also estimate the total number of ionizing pho- 
tons per baryon required to maintain the ionized IGM 
between Zcnd and z, 

- xMz')dt_^^,^ (29) 

^rcc 



where Zcnd = 6 is the redshift of the end of the reioniza- 
tion, and if we assume the helium is also totally singly 
ionized at Zend = 6 we can get iVcnd — 1- In fig. [31 we 
show the N[^^^{z) at different z for three /* cases. We 
find the number keeps going up until around z = 15 and 
becomes constant ~ 2.5 at higher redshift for three cases 
of the reionization histories. 

5. THE NEAR-IR EBL INTENSITY FROM REIONIZATION 

The mean cosmic infrared background can be esti- 
mated by 



dz 



V{z)jy{z) 



H{z) (1-t-z) 



(30) 



where v(z) = (1 + z)vohs, and we take Zmin — 6 and 
•^max = 30 in the calculation. This redshift range can 
fully take account of the emission from the Pop III and 
early Pop II stars, which is rcdshifted into the near-IR 
band. We take this to be the form of 



3u{z) = /pjp°p"i(z) + (1 - h)ru°'"\z), 



PopII/ 



(31) 



where and are the comoving specific emis- 

sion coefficients 



(32) 



where l\, is the luminosity mass density at v, {t'D is the 
mean stellar lifetime of each of the stellar type and tp{z) 
is the comoving SFRD given by Eq. (125]) . 

In Fig.m we show the spectrum of the near-IR back- 
ground light intensity from both Pop II and Pop III stars. 
We assume that the reionization is ending around z = 6 
and integrate up to z = 30 to determine I^. Similar to 
Fig. [H we plot the total spectrum from the stellar neb- 
ula and the IGM for both Pop II and Pop III stars. The 
contributions from different sources we consider are also 
shown in colored thin lines. 

Here we still take the same value for and fcsc as in 
Fig. 1 related to the li, calculation. Similar to the lu- 
minosity mass density, the spectrum from stellar nebula 
is much larger than that from the IGM, and the stel- 
lar spectrum is higher for Pop II stars while the "back- 
ground" spectrum is higher for Pop III stars. However, 
different from the luminosity mass density, we now find 
the spectrum from Pop II stars is larger than that from 
Pop HI stars. This is basically because the typical life- 
time of Pop II stars is longer than that of the Pop HI 
stars. 

6. ANGULAR POWER SPECTRUM 

The angular cross power spectrum of the infrared emis- 
sion at observed frequencies v and i^' for a multipole £ 
is 



dzl^ 

dz 



-j j<^iz)ji^'iz)Pgf,{k,z), 

(33) 

where x is the comoving angular diameter distance, 
a=(l + z)^^is the scale factor, and jv{z) is the mean 
emission per comoving volume at frequency v and red- 
shift z. If we just take account of the flux lower than a 
upper cut S'cut, and then ju{z) can be written as 



j,{z) ^(l + z) 



dS S 



d^N 
dSdz 



(34) 



Here S is the source flux and N is the number of sources. 
This quantity is just the comoving specific emission co- 
efficient we derive in the last section. The Pgg{k,z) is 
the galaxy power spectrum at wavenumber k — £/x and 
redshift z, and we will make use of the model of halo 
occupation distribution to calculate the Pgg(k,z). 

6.1. First Galaxy Clustering 

To calculate the Pgg we extend the linear theory model 
of Cooray et al. (2004) and make use of the halo occu- 
pation distribution (HOD) for first-light galaxies during 
reionization. The galaxy power spectrum can be written 
as 



Pgg(fc, z) = P^^{k, z) + P^^ik, z), (35) 

Where Pgg and Pg^ denote the power spectrum con- 
tributed by galaxies in a single dark matter halo and 
galaxies in two differ ent dark ma tter halos respectively. 
Then we can write (C oorav fc: Sh cth 2002) 



P 



Ih _ 



dM 



dn (A^gal(iVgal- 1)) 

dM nl^. 



uP{M,k), (36) 



P 



2h _ 



-P 



lin 



dM b{M,z)4^^-^u{M,k) 

^ ' 'dM Ugal ^ . 



.(37) 
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Fig. 4. — The spectrum of the near-IR EBL intensity vs. the observed frame wavelength Ao^s fof the Pop II and Pop III stars integrated 
from z = 6 to 30. We take /» = 0.03 and /esc = 0.5 for both Pop II and Pop III cases. We find the spectrum of the stellar nebula is 
much larger than that from the IGM for both Pop II and Pop III stars which is similar with what is shown in Fig. [T] However, unlike the 
luminosity density case, the spectrum from the Pop II stars is larger than that from the Pop III stars. This is basically because the typical 
lifetime of the Pop II stars is larger than that of Pop III stars. 



Here M is the halo mass, dn/dM{M,z) is the halo 
mass function, u{M, k) is t he Fourier transform of the 
NFW halo density profile (jNavarro et all I1995D . p = 
1 wh en (N^p,](N^p,] — 1)) < 1 and p = 2 other- 
wise (iCoorav fc ShethI 120021 ) . b{M,z) is the halo bias 
(jSheth fc Tormenlll999D. and P\i„(k, z) is th e linear mat- 
ter power spectrum (jEisenstein fc Hulll997D . The n^ai is 



the mean number density of galaxies, which is given by 



(38) 



The (A^gai) is the mean number of galaxies in a halo with 
mass M, which is the sum of number of central galaxies 
and satellite galaxies (jZheng et al.ll2005() 



(iVgal) = (A^con) + (iV.at), 



where we define 



1 + erf 



logAf - logM„ 



(39) 



(40) 



and 



1 + erf 



\ogM - logMo 

Clog A/ 



/ M 



(41) 

In this definition, the Mmin denotes the mass of a halo 
that has 50% probability to host a central galaxy, and 
ciogM is the transition width. For the satellite galaxies, 
Mq is the truncation mass for satellites, Mgat is the nor- 
malization mass and as denotes the slope of the power- 
law relation about the halo mass M. We assume Mq is 
always larger than M„iin since there should not be satel- 
lites without central galaxy, and assume Mq = 2Mmin. 
We take Afgat = 15Mi„in, criogM = 0.3, and as = 1-5 in 
this paper. If assuming a Poisson distribution for satel- 



lite galaxies, we can get 

(A^gal(iVgal-l)) ~2(7Vsat)(^c, 



(42) 



This expression could take account of the case < 
(Nccn) < 1 and is consistent with our definitions for the 
(A^ccn) and (A^sat). 

6.2. Poisson Fluctuations 

The clustering measurements are affected by the Pois- 
son fluctuations associated with the shot-noise caused 
by the discrete and finite number of galaxies from which 
clustering is measured. Assuming a Poisson distribution 
the ^-independent shot-noise power spectrum is 



Cshot 
e — 



Sou 







(43) 



To estimate Cf^°^, we first define the luminosity mass 
density for the mass of the dark matter halos at frequency 

V 



r 



where = L^/IVL^ is the luminosity mass density for 
the stellar mass discussed in Section 3. We derive the 
3-D shot-noise power spectrum by assuming is pro- 
portional to the halo mass M, i.e. is independent on 
M 



P: 



;hot 



I'' 



(44) 



Then the 2-D shot-noise power spectrum can be written 
as 



shot 



Pf°'(z). (45) 
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6.3. Band-Averaged Intensity Power Spectrum 

For a specific near-IR observation with a band fre- 
quency from vi to 1^2, we can define a band-averaged 
luminosity mass density as 

Z = — / dvl,, (46) 

Jui{l+z) 

wfiere Ai/ — V2 — v\ is the bandwidth. Tlien we can 
derive the band-averaged comoving specific emission co- 
efficient j{z) using Eq. ([5^ and the 3-D shot-noise power 
spectrum P'^hot ^gjj^g pi|) respectively. 

Finally we find the band-averaged angular cross power 
spectrum and shot-noise power spectrum are 

C, = dz j{z)j{z)P,,{k, z), (47) 

and 

respectively. Note that here we have a factor instead 

of o? in C^" and C^^ an d this dependence has b een 
explained in the Appendix of iFernandez et all (j201Cl( ). 

7. RESULTS AND DISCUSSION 

In this Section, we first estimate the infrared back- 
ground intensity and then discuss the angular power 
spectrum as derived previously. We also compare our 
estimation with the observational data and discuss the 
dependence of the result on the parameters in the model. 

7.1. Near-IR EBL from Reionization 

In Fig. [51 we show the spectrum of the near-IR inten- 
sity vs. the observed frequency for three cases of the 
reionization history with r = 0.077, 0.090 and 0.099 
corresponding to /, — 0.02, 0.03 and 0.04. Here the 
spectrum is the total spectrum of the sum of that from 
both stellar nebula and the IGM, which are calculated 
by putting Eq. 1^ into Eq. Also, we set /esc = 0.5 
for both Pop II and Pop III stars. The right panel of 
Fig. [5] shows three cases with esc varied at fixed /, . In 
both panels, for comparison, we plot the observational 
data in terms of the excess of EBL relative to the inte- 
grated galaxy light (IGL) of known galaxy populations 
at low redshifts (Madau & Pozzetti 2000). The excess 
EBL data plotted here are the same as those shown in 
Santos et al. (2002) and involves measurements mainly 
from DIRBE with various models for zodiacal light and 
Galactic stellar contribution. 

We find the emission from the Pop II stars dominates 
the spectrum for all of the three cases. As we have just 
discussed in the last Section, this is because the lifetime 
of the Pop II is longer than the Pop III stars. In the spec- 
trum, the shorter and longer wavelength parts are mainly 
contributed by the "background" spectrum from Pop III 
stars while the medium part by Pop II stellar spectrum, 
so if just the Pop III stars get longer lifetime, only the 
"background" spectrum can be effectively reinforced in 
the total spectrum. 



In any case, regardless of assumptions related to the 
stellar type, we find that the EBL intensity from reioniza- 
tion is no more than 0.4 nW m~^ sr~^. Such an intensity 
is significantly smaller than the previous predictions that 
attempted to explain almost all or a significant fraction 
of the excess EBL seen in DIRBE data relative to IGL 
estimates. An intensity larger than about 2.5 nW m^^ 
sr~^ in the J-band could be in conflict with metal produc- 
tion considerations and the X-ray background (Madau & 
Silk 2005), though they do not necessarily require high 
efficiency factors to generate the required star formation 
(Fernandez & Komatsu 2006). The difference between 
our calculation and the previous ones is that we primarily 
require the reionization model to generate a reionization 
history consistent with the WMAP optical depth. This 
limits the number of H-ionizing photons per baryon dur- 
ing reionization to be less than 3. Previous estimates ig- 
nored such a constraint and either focused on explaining 
all of the DIRBE excess (Santos et al. 2002; Salvaterra 
& Ferrara 2003) or amplitude of the measured near-IR 
background anisotropics (e.g., Fernandez et al. 2010). 

In Fig. [Bl we show the total near-IR background inten- 
sity spectrum for different transition redshifts zt from 
Pop II to Pop III stars with increasing redshift. Here 
we take four transition redshifts Zt = 10, 15, 20 and 30 
using Eq. (|26l) . To maintain the same optical depth to 
electron scattering r — 0.090, we find a higher intensity 
is predicted for higher z^. Correspondingly should be 
increased to values of 0.040, 0.043 and 0.044 for zt = 15, 
20 and 30, compared to — 0.03 for Zt — 10. If Zt = 30 
(red dotted line) the integrated intensity is about ^ 3 
to 4 times greater than the case for zt = 10 (blue solid 
line) when Aobs > 2/im. Also, the peak of the integrated 
intensity spectrum moves to longer wavelengths as zt is 
increased. This is because the Pop HI stars are generally 
hotter than the Pop II stars, which can produce more 
ionizing photons. When the transition redshift is higher, 
the longer the Pop II stars dominate the Universe, and 
less ionizing photons are produced when compared to the 
case for a low zt. To keep the reionization history un- 
changed, we need more Pop II stars to generate enough 
ionizing photons. This results in a near-IR background 
intensity that is higher. 

7.2. Bright-end Galaxy Luminosity Functions 

In order to relate our galaxy population responsible 
for reionization to the observations, we also compare our 
model to observations of z > 6 galaxies, focusing on the 
UV LFs. To obtain the UV LF, instead of the occupation 
number which is luminosity independent, we make use 
of the conditiona. 1 luminosity function (CLE) approach 
(I Ya ng et al.l[2003L iCoorav fc Milosavlievidl2005t iCooravl 
2005). To compute the CLE we map galaxy rest- frame 
UV luminosity to halo mass with some scatter added 
similar to the case of low-redshift galaxy populations 
(Cooray 2005). The total luminosity of a halo is taken 
to be 

Ltot(M,z) =/,(z)/,^Af, (49) 

and we assume that this total luminosity can be ascribed 
to the central galaxy when Mmin < M < Mgat , following 
the' earlier description related to the occupation num- 
ber. When M > Mgat we introduce satellites with cen- 
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Fig. 5. — The total near-IR background intensity spectrum for three assumed values of /, at fixed /esc (left panel) and the same with 
/esc varied at fixed /* (right panel). In the left panel, the dashed, solid and dotted lines denote the three cases of reionization history with 
T = 0.077, 0.090 and 0.099. The models in the right panel also range in t between 0.07 and 0.1. The data points with error bars are from 
Santos et al. (2002) and show the "DIRBE excess" defined as the difference between total DIRBE background and the integrated galaxy 
light at each of the wavelengths. The points shown with arrows are strictly upper limits since the low-end of the error is consistent with 
zero for them with significance less than 2a. 
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Fig. 6. — The total near-IR background intensity spectrum for 
transition redshifts zt of Pop II and Pop III stars between 10 and 
30. To maintain the same optical depth to electron scattering of 
T = 0.090, a higher intensity is predicted as zt is increased. For 
example the integrated intensity when zt = 30 (red dotted line) 
is ~ 3 to 4 times greater than when zt = 10 (blue solid line) for 
Aobs > 2^111. 

tral galaxy luminosity kept fixed at Ltot(-^^sat, -z)- How- 
ever, when comparing to the existing measurements we 
found that all of the rest-UV LF measurements are in the 
range where central galaxies dominate the LF and thus 
our comparison to the measured LFs is independent of 
assumptions related to the exact form of the satellite oc- 
cupation number or conditional luminosity function. 
To compare with existing rest-UV LF measurements, 



we convert the luminosity of each galaxy to the AB ab- 
solute magnitude via the relation Mab = — 2.51ogioiiy + 
5.48 0. In the left panel of Fig. [71 we show the 
SFRD as a function of redshift derived from Eq. (f25|) 
for three reionization histories with t = 0.077, 0.090 
and 0.099 which are obtained by setting /* — 0.02, 
0.03 and 0.04, respectively. The red points are 
the data from HUDFq 9-HERS-hCANDELS observations 
(jBouwens et al.ll2011bl ). We find the SFRD of the three 
cases are higher than the existing measurements, espe- 
cially at the high redshifts. This difference is mainly due 
to the fact that existing SFRD estimates are limited to 
galaxies with Muv < —17, while the bulk of the reion- 
ization UV density budget is contained in the galaxies 
at the faint-end of the LF. This is especially the case at 
high redshift since the faint-end slope of the LF is steep 
with values reaching close to -2 already. 

In the right panel of Fig. [71 we show our rest-UV LF 
(corresponding to Arcs = 1600 A) at z =6, 7 and 8. 
The central thick solid curves are the luminosity func- 
tion derived from our default model with — 0.03 and 
T = 0.090, and the thin dotted lines are obtained by 
/* = 0.02 (lower) and 0.04 (upp er) which could matc h 
the la errors of the data given in (jBouwens et al.l[2011a[ ). 
The three values of the faint-end slope a = —2.5, -2.0 and 
-1.5 are shown, which indicates the slope of our model is 
between -2.5 and -2.0. Also, we find the star formation 
time scale tsF is around 6 x 10^ yrs at ^ = 6, 7 and 8 
when we calculate the luminosity mass function here. 

We also explore the dependence of the galaxy bias on 
the Mab in Fig. [H We can define a effective galaxy bias 
here from the HOD model as 

5eff (z) = r^^'^ dM-^biM, z)^^ , . (50) 

http://www.ucolick.org/ cnaw/sun.html 
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1600.AB 

Fig. 7. — Left: The SFRD as a function of the redshift derived from Ea. H25| l for three reionization histories t = 0.077 (blue dashed Une), 
0.090 (blue solid line) and 0.099 (blue dotted line) which are obtai ned by setting ft = 0.0 2, 0.03 and 0.04, respectively. The red points are 
the measurements from HUDF09+ERS+CANDELS observations l IBouwens et al.ll2011bl ). Right: The estimated UV luminosity function 
at 2 = 6, 7 and 8 with Arcs = 1600 A. The central thick solid lines are derived from our model with /, = 0.0 3 at z = 6, 7 and 8, an d the 
thin dotted lines are for /» = 0.04 (upper) and 0.02 (lower) which can match the la error of the data given in lIBouwens et al.l l2011a!). The 
faint-end slope a at three values are also shown in black solid lines for comparison. 

tering power spectra with non-linear power spectrum 
from the HOD model are in solid lines. The left panel 
shows the case with M^in, the minimum mass to host 
a galaxy, is changed from 10^ to 10^ Mq. For com- 
parison, we also plot the linear power spectrum which 
has a turnover around £ = 10^. When calculating the 
clustering power spectrum, we set the parameters of 
the HOD with Mmin values as listed in the figure with 
Msat = ISMmin- The corresponding values on the optical 
depth to electron scattering are also listed in the figure. 
We note that the shot-noise amplitude is larger for the 
case with Mmin = 10^ Mq in comparison to the case with, 
say, Mynin = 10^ M©, though in those two cases the clus- 
tering amplitude is higher with Mmin — 10^ M©. This is 
because the shot-noise amplitude is sensitive to the sec- 
ond fiux-moment of the number counts. By keeping the 
minimum mass higher we force the overall counts to be 
restricted to brighter sources than the case with a lower 
minimum halo mass. On the other hand the clustering 
power spectrum refiects the total background intensity. 
With the minimum mass lowered, both the overall num- 
ber density of galaxies and the background intensity are 
increased. 

The right panel of Fig. |9] shows the case where we vary 
/esc to highlight the fact C; amplitude is inversely pro- 
portional to /ese- However, one cannot arbitrarily reduce 
/esc to a small value since this results in a low optical 
depth to electron scattering. If M^nin ~ 10® to 10^, such 
that one does not need to be concerned of effects due to 
feedback negatively impacting the formation of galaxies 
in lower mass halos, then we find that /esc must be at 
the level of 0.3 or more, under the assumptions related 
to the IMFs for PopH and PopIII stars we use in this 
paper. 

In Fig. [TOl we plot the Ci at different redshift bins. As 
can be seen, the main contribution of the Cf comes from 




Kb 

Fig. 8. — The effective galaxy bias vs. AB magnitude A/ab 
2 = 6, 7 and 8. We find the bias increases at higher redshifts and 
lower Mab- 

In the plot we show the bee as a function of Mab at 

z =6, 7 and 8 for the Mab = -22 17. We find that 

the galaxy bias increases as the redshift increases and 
decreases as the Mab increases. The reason is obvious 
that the galaxy number density figai defined by Eq. ([38]) 
becomes smaller at higher redshift and bigger at larger 
Mab (the larger Mab means smaller halo mass M). 

7.3. Anisotropy Power Spectrum 

In Fig. |9]we show the near-IR background anisotropy 
angular power spectrum at Aobs = 3.6 /Ltm. The clus- 
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Fig. 9. — The near-IR background intensity anisotropy power spectrum at Aq 



: 3.6 fim. The solid lines are the total clustering power 



spectrum with the non-linear power spectrum from the HOD model. Left: We show a variety of predictions with Mmin related to the 
occupation number taken to be values shown in the plot. In all of these cases we take Msat = ISMmin- As Mmin corresponds to the 
minimum halo mass to host a galaxy, the reionization history is also changed and we list the optical depth to electron scattering for the 
four cases we have studied here. For comparison, the linear power spectrum is shown in a dashed line for one of the cases, which have a 
turnover around £ = 10"^ (Cooray et al. 2004). The shot-noise power spectra are in dotted lines. Note that the shot- noise level is higher 
when Afnjin is larger. Right: Clustering predictions as a function of /esc; Note that the amplitude of the clustering power spectrum is 
higher with a smaller value for /esc- However, the optical depth to electron scattering is lower with a small value of /esc. 

the lowest redshift range of 6 to 10. For all practical 
purposes one can assume that the near-IR background is 
probing the end of reionization not the first objects to 
form in the universe at the beginning of reionization. 




Fig. 10. — The angular power spectrum at 3.6 fim for different 
redshift bins. The power spectrum dominates over the redshift 
range of 6 to 10. The model shown here is the default case with 
r = 0.090 with /esc = 0.5 and /, = 0.03. 



To compare with the total near-IR intensity spectrum, 
we plot the ratio of the square root of Ci to the vli, in 
Fig. [11] We integrate from z = 6 to 30 with the same 




Fig. 11. — The ratio of rms fluctuations to the total intensity at 
different wavelengths X^hs = 1-6, 2.4, 3.6 and 4.5 fim. We find the 
near-IR background fluctuations are around 10% compared to the 
background intensity from £ = 10"^ to 10*. 



HOD model described above. For easy comparison, we 
restrict the predictions to be with Mmin — 10^ Mq such 
that T ~ 0.090 to remain consistent with WMAP 7-year 

result. We find the ratio of y^£{i + 1)Ci/2tt to vly is 
about 10% for 10^ < € < lO'' and the clustering fluc- 
tuations amplitude is below 100% of the intensity at all 
angular scales. This shows that large fluctuations of the 
background intensity is not expected and the background 
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behaves in a manner that is smooth and not clumpy as 
in the case if spatial variations are dominated by rare, 
bright galaxies. 

In Fig. [121 we show the near-IR background inten- 
sity power spectrum (solid curves), shot-noise power 
spectrum (dotted curves) with a comparison to exist- 
ing observational measurements at several wavelengths : 
Aobs = 1.6, 2.4, 3 6, and 4.5 urn (iThompson et"aIl l2Q07l: 
iMatsumoto et all 120111: iKashlinskv et al.l l2012(l Note 
that for the data from Kashlinsky et al. (2012) at 3.6 
and 4.5 /im, we removed the IRAC beam transfer func- 
tion Bi so the comparison to theoretical predictions is 
easier. An accounting of the beam is essential since real- 
istic instruments do not have perfect resolution causing 
a loss of power on the smallest scales. This effect can be 
explicitly measured from the point spread function. We 
calculate the beam transfer function Bi by taking the 
power spectrum of the measured point spread function 
for the Spitzer/IRAC instrument, which is publicly avail- 
able and compute Ci — >■ Ci/ Bi (for more information see 
Smidt et al. 2012). 

We also show the results from the simulation in 
[Fernandez et al.l (|2012[ ) for comparison. The red solid 
and magenta dashed lines are for the simulation with 
mass resolution M^i^ = 10* and Mmin = 10^ M©, re- 
spectively. These two cases are what is described in their 
work as not having a suppression of the small mass halos 
(ie star-formation present in halos with mass between 10* 
and 10^ Mq, where one expects photo- ionization heating 
to suppress star-formation) . The reionization histories of 
these two cases involve PopII stars with /esc = 0.1. 

For more clarity, we show the square root of the Ci 
at ^ = 3 X 10'^ and ^ = 10'* together with the data as a 
function of the wavelength in Fig. [131 The data points 
are from the data sets shown in Fig |12l which are around 
^ = 3 X 10'^ or ^ = 10*. We find the curve of the square 
root of the Ci has similar shape as Vohslu, but the am- 
plitude of rms fluctuations in our models is lower than 
the existing measurements. We capture the uncertain- 
ties related to /esc, and M,„i,i by considering a low 
and high range for our prediction related to Ci with r 
falling within the 1 a uncertainty range of the WMAP 7- 
year result when these parameters are varied. Even with 
parameter uncertainties accounted for, we find that the 
existing measurements are at least an order of magnitude 
larger than our model predictions. 

We attempted additional model variations but failed to 
find a scenario where r is consistent with WMAP 7-year 
result and the existing LFs leading to a model consistent 
with existing near-IR background. One can, in princi- 
ple, model fit the near-IR fluctuation power spectra by 
increasing the photon output of first galaxies. This re- 
sults in an optical depth that is higher than the WMAP 
value and a LF that has brighter galaxies than observed 
in existing deep HST/WFC3 images. The existing mea- 
surements require a background intensity that is at least 
3 nW m^^ sr~^ at 1.6 /^m so that the model predictions 
become consistent with measurements within the la un- 
certainties of the measurements. Since 0.3 nW m~^ sr~* 
at 1.6 /im is associated with 2.5 H-ionizing photons per 
baryon during reionization, if reionization is to explain 
the background fluctuations then we are dealing with a 
situation where ^ 25 H-ionizing photons per baryon are 



present. A possibility is to introduce a spectrum for the 
emission that has a rest-frame cut-off in UV at wave- 
lengths shortward of the Lyman-limit that is not asso- 
ciated with reionization but by the emission mechanism 
itself. Unfortunately we have not been able to come up 
with such an emission spectrum. 

As discussed in Helgason et al. (2012) faint, low- 
redshift galaxies do not have a clustering shape consis- 
tent with the anisotropy power spectrum measurements. 
The existing measurements show a clear excess in clus- 
tering at 30 arcsec to few tens arcminute angular scales 
that differ from faint galaxy clustering power spectrum. 
This difference is statistically significant and one can only 
explain at most 20% of the near-IR fluctuations to be 
associated with faint galaxies, reducing the previous es- 
timate in Cooray et al. (2007; Chary et al. 2008) that 
suggested a contribution at the 50% level. We are now 
left with an unexplained set of measurements since nei- 
ther the low redshift faint galaxies nor the high redshift 
reionization galaxies can explain them. One possibility is 
that the existing fluctuation excess in the data has a non- 
astrophysical origin, perhaps either involving systematics 
in the data or fluctuations in the Zodiacal light. While 
arguments have been made that these effects are insignif- 
icant with a set of well-coordinated multi-wavelength 
measurements we plan to further address the origin of 
near-IR fluctuations in our upcoming papers. 

8. SUMMARY 

The UV emission from stars which are formed in the 
early Universe from 2; = 6 to 30 can contribute to the 
near-IR background light. By measuring the intensity 
and anisotropies of the near-IR background, we can in- 
vestigate the properties of these early stars and the epoch 
of reionization. In this work we discuss several sources 
that contribute to the near-IR background intensity, in- 
cluding the emission of stars, Lya, free-free, free-bound 
and two-photon. We first estimate the frequency spec- 
trum separately for Pop II and Pop III star with the 
redshift range z = 6 ^ 30. By using the initial stellar 
mass spectrum and the fitting model of time-averaged 
hydrogen photoionization rate we calculate the luminos- 
ity mass density. We find that although the luminosity 
mass density of the Pop III stars is a bit larger than that 
of Pop II stars, the near-IR intensity spectrum from the 
Pop II stars is stronger than that from Pop III stars, 
which is caused by the longer lifetime of Pop II stars. 

In order to check the consistence of our stellar model 
and the reionization history, we derive the hydrogen 
reionization fraction at different redshifts and calculate 
the optical depth by assuming the Pop II stars are dom- 
inant for z < 10 while Pop III stars for z > 10. We find 
that if we set the star formation efficiency /* = 0.03, the 
Universe would be totally reionized around z = 9 with 
the optical depth r — 0.090, which is well consistent with 
the result from the WMAP 7- year data. Also, we explore 
the other possible models with /» — 0.02 and 0.04 and 
get r = 0.077 and 0.099, respectively. The total number 
of the ionizing photon per baryon required to maintain 
the ionized IGM, -/V^n(z), is also evaluated for the three 
cases, and we find the NP^^^z) becomes a constant with 
a value around 2.5 after z ~ 15. 

To compare with existing bright-end luminosity func- 
tion measurements, we evaluate the UV luminosity tunc- 
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Fig. 12. — The comparison of our Ci with the current measurements at different wavelengths. The clustering power spectrum are in blue 



solid (/* = 0.03 and /esc = 0.5), dashed (/, = 0.02 and /esc = 0.9) and dotted (/* = 0.04 and /esc = 0.1) lines with M„ 



10*^ Mq in 



our HOD model. The shot-noise power spectrum are in blue dash-dotted lines. The observational data of shot-noise are shown in black 
dash-dotted lines. For the data of Kashlinsky et al. (2012) at 3.6 and 4.5 lira, we r emov ed the beam function so that a direct comparison 
can be made to model predictions. The simulation results from [Fernandez et al.l II2012I ) are also shown with minimum halo mass set at 
Mjnin = 10* (lower red solid) and Mj^in = 10^ Mq (upper magenta dashed), without a suppression related to the star- formation in halos 
between the range of 10* and 10^ Mq where photoionization effects become important. The shaped region in each of the two lower panels 
is an estimate of the residual anisotropy power spectrum signal expected from low-redshift faint galaxies following Helgason et al. (2012). 
These two curves clearly show the presence of excess clustering at ^ ~ 10'^. While this excess has been suggested as originating from 
reionization (Kashlinsky et al. 2012), we do not find this to be the case as our predictions are lower. 



tion from our model at Arcs = 1600 A for z —Q, 7 and 
8 from Mab — ~22 to —17. We compare to the mea- 
surements of Bouwens et al. (2012). We find our derived 
luminosity function is consistent with the data with the 
slope of the faint-end a ~ —2. This is a steep slope 
and existing measurements do suggest that the slope is 
steeper than for LFs at low redshifts. We then define 
an effective galaxy bias h^s from the HOD model, and 
find the 6eff becomes bigger at high redshift and bright 
luminosities because of the lower number density of the 
galaxies. 

Finally, we calculate the angular power spectrum of 
the near-IR background by making use of a halo model. 
The non-linearities provided by the 1-halo term increase 



the clustering strength at multipoles greater than about 
^ ~ 2 X 10'^. The suggest turn-over at this scale with the 
linear power spectrum alone in Cooray et al. (2004) is no 
longer present. We find the shot-noise power spectrum 
of the Pop III stars is greater than that of Pop II stars 
because of the larger luminosity mass density for Pop III 
stars. By making use of the stellar population evolution 
model, we calculate the near-IR anisotropy power spectra 
at different wavelengths and compare to the observa- 
tional data. We find our results are lower than the obser- 
vational data by at least an order of magnitude. There 
are now strong arguments that the near-IR anisotropics 
cannot originate from low redshift faint galaxies (Helga- 
son et al. 2012). We have failed to explain the alterna- 
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Fig. 13. — The square root of the Ce at £ = 3 X 10'' and £ = 10* compared to the observational data at different wavelengths Aq^s = 1-6, 
2.4, 3.6, and 4.5 fim. The blue solid, dashed and dotted line are the square root of the Cg with /, = 0.03 and fcs c = 0.5, /« = 0.02 and 
/osc = 0.9, and /* = 0.04 and fcsc = 0.1, respectively. The data points are from the data sets shown in Fig. 1121 around £ = 3 X 10"^ or 
£ = 10"'. 



tive origin of near-IR background anisotropics involving 
galaxies present during reionization, contrary to sugges- 
tions in the literature (Kashlinsky et al. 2012). In future 
works, using additional measurements with Spitzer and 
Hubble/WFC3, we plan to further discuss the near-IR 
fluctuations and to explain if the origin is astrophysical 
or whether it is associated with yet-unknown systematic 



effect in the data. 
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